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Abstract. We develop a new model-independent method to probe the constancy of the 
speed of light c. In our method, the degeneracy between the cosmic curvature and the speed 
of light can be eliminated, which makes the test more natural and general. Combining the 
independent observations of Hubble parameter H(z) and luminosity distance dz(z), we use 
the model-independent smoothing technique, Gaussian processes, to reconstruct them and 
then detect variation of the speed of light. We find no signal of deviation from the present 
value of the speed of light cg. Moreover, to demonstrate the improvement in probing the 
constancy of the speed of light from future experiments, we produce a series of simulated 
data. The Dark Energy Survey will be able to detect Ac/co ~ 1% at ~ 1.50 confidence 
level and Ac/co ~ 2% at ~ 30 confidence level. If the errors are reduced to one-tenth of the 
expected DES ones, it can detect a Ac/co ~ 0.1% variation at ~ 20 confidence level. 
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1 Introduction 


The constancy of the speed of light is one of the most fundamental and recognized physical 
properties and plays a crucial role in basic physical laws such as Maxwell equations, special 
and general relativity and many others. Although many measurements of the speed of light 
have been performed since 1675 and ended with very accurate measurements today, some 
people argue that the speed of light c may be dynamical and has been varying in the past. 
Such an idea has attracted a lot of interest recently because it can provide an alternative 
solution of horizon and the flatness problems in the non-inflationary cosmology. The the- 
ories of varying speed of light (VSL) have been considered [1-4]. However, a comparison 
of those theories with experimental data seems to be still missing. Recently, V. Salzano et 
al. introduced a model-independent method to measure the speed of light through Baryon 
Acoustic Oscillations (BAO), and then it was followed by some newly updated results about 
its application to the forecast data [5]. In their paper, based on a spatially flat Friedmann- 
Robertson-Walker (FRW) universe they derived a relation among the speed of light, the 
angular diameter distance and the Hubble function, 


Da(zm)- H(zm) = c(zm), (1.1) 


where zm is the redshift at which the angular diameter D4(zm) is maximal. They used the 
Gaussian Processes (GP) to reconstruct D,4(z) and H(z) and then found the redshift zm. 
Finally, they used Eq. (1.1) to obtain the speed of light at zm and probed the constancy of 
the speed of light. However, since zm is not covered by observational data sets, they only 
used mock data to test the accuracy of the method. More recently, they showed that Square 
Kilometer Array (SKA) can detect a 1% variation in the speed of light at 30 level [6]. But 
smaller signals are hardly detected by already-planned future galaxy surveys. 

It seems that the method has two drawbacks. The first is that the speed of light can be 
measured only at one redshift zm, which is extremely limited by the redshift range covered 
by the data. The second and most important is that they ignore the cosmic curvature which 
is degenerated with the speed of light. Although the value of the cosmic curvature given 
by current data is very small, if we want to probe even smaller variation of the speed of 
light, it may be problematic. As a model-independent method, it is more general to include 
the cosmic curvature. In this paper, we develop a new and more general model-independent 
method which can dodge the cosmic curvature. We use the luminosity distance dz instead 
of the angular diameter distance D4 to test the speed of light, so that we can use the large 
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sample of real data such as cosmic chronometers (CC), BAO and SNela Union 2.1 or joint 
light-curve analysis (JLA) to probe the constancy of the speed of light over a wide range of 
redshift. We will also produce a series of mock data to test the accuracy of our method. 
This paper is organized as follows. In Sec. 2, we introduce the theoretical method to 
probe the constancy of the speed of light. In Sec. 3, we apply GP method to the test using 
two independent data sets: CC+BAO and SNela Union 2.1 (and JLA). Furthermore it is 
followed by a series of simulated data tests. We give discussions and conclusions in Sec. 4. 


2 Method 


For a FRW universe, given that the speed of light is a function of time, the line-element can 
be expressed as 


dr? 


2 2 2 2 
ds* = —c (t)dt +a (t) I Kr 


+ r?(d6? + sin?bde’) | , (2.1) 
where K = +1,—1,0 corresponds to a closed, open and flat universe, respectively. Note 
that we write the metric with a preferred proper comoving time [7, 8]. We define D(z) = 
(Ho/co)(1 + z)~'dz(z) as the normalized comoving distance, where co is the value of the 
speed of light today, and dz(z) is the luminosity distance. Given that the speed of light is 
time dependent, we should reevaluate the luminosity distance from the definition: 


L 
dœ = a 
LAr’ 


(2.2) 


where L, is the absolute luminosity of a source and F is an observed flux. Setting r = sin x 
(K =+4+1), r=y (K = 0), and r = sinh y (K = —1) in Eq. (2.1), the 3-dimensional space 
line-element is expressed as 


do? = dy? + (fr (X))? (db? + sin26de?), (2.3) 
where 
siny (K = +1), 
fro) = 4 X (K = 0), (2.4) 
sinhy (K =-1). 
The function (2.4) can be written in a unified way: 


1 
V/-K 


We can easily check that x is the comoving distance: 


_ 1 f cf) z 
mh, EO” me 


Note that the observed luminosity Lo (detected at x = 0 and z = 0) is different from the 
absolute luminosity L, of the source (emitted at the comoving distance x with the redshift z). 
The flux F is defined by F = Lo/S, where S = 4r ( fg (x)}? is the area of a sphere at z = 0. So 
the luminosity distance in Eq. (2.2) yields d2 = (fx(x))?Ls/Lo. If we write the energy of light 
emitted at the time-interval At; to be AE}, then the absolute luminosity Ls = AF, /At,, 


f(x) = sinh(V—K x). (2.5) 
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and similarly the observed luminosity is given by Lo = AEp/Atg. Note that here the speed 
of light is time dependent, so AF, /AEp = ™/v = (c1/A1)/(co/Ao) = E(z)(1 + z), where 
êlz) = c(z)/co. Moreover, Ato/At; = vı /v = G(z)(1 + z). Hence we find 


dr = e(z) fx(x)(1 + 2). (2.7) 


Then combine the definition of D(z) and Eqs. (2.5)-(2.7), we can derive 


Fe ee eer ore / . Ha, (2.8) 


WOK 


where E(z) = H(z)/Ho, and Qe = —Kco?/Hpo? is the dimensionless curvature density 
parameter at z = 0. Differentiating Eq. (2.8) with the redshfit z, we can get 


o, — LOD = ¢@POPE - &) 
1 a . 


D(A (2) a 


From Eq. (2.9) we can see that there is degeneracy between VSL and curvature which is also 
pointed out in Ref. [6]: the possible detection of a signal might be equally interpreted as 
“VSL + null curvature” or “constant c(z) + curvature”. Thus, to eliminate Qg, using the 
second derivative of D(z) we obtain 


O= @(z) 4+ A(z)é(z) — B(z)é(z)é (z) 


+ M(z)e?(z) — N(z)é(z)é"(z), (2.10) 
where 
A(z) = [D"(z)E(z) + D'(z)E"(z)| D(z)E(2) 
— D(z) E* (2), (2.11) 
B(z) = E? (2) D(2)D' (2) + D”(2)E(2)E' (2), (2.12) 
and 
M(z) = 2N (z) = 2D? (2) E? (2). (2.13) 


Eq. (2.10) is valid for any form of c(z) which is the function of time no matter what the 
value of the cosmic curvature is. If we assume the constancy of the speed of light, thus 
êlz) = c(z)/co = 1, then Eq. (2.10) becomes 


T(z) =14+ [D"(2)E(2) + D'(2)E'(2)| D(2Z)E(2) 
— D?(z)E*(z) =0, (2.14) 


it always holds if the speed of light is constant. Any deviation of T from 0 at redshift z, 
will indicate c(z,) is different from co. Thus we can test the speed of light at every redshift 
we want. This will make our test more flexible so that our method will not be limited to 
probing the constancy of the speed of light only at redshift zq as obtained in [5, 6]. On the 
contrary, we can choose the redshift where the quality of the data is better so that we can 
improve the precision of the test for the mock data. 

In summary, we propose a null test Eq. (2.14) to probe the constancy of the speed of 
light in a more general way, which is cosmological model-independent. The most important 
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point is that we can dodge the cosmic curvature. Thus we need not assume the value of 
the cosmic curvature. It is emphasized that the only hypothesis is the FRW metric of the 
background. In addition, we would like to mention here that once reconstruct A(z), B(z), 
and M(z), we can obtain the dependence of the speed of light on the redshift through solving 
(2.10). However, in this paper we pay attention on testing the constancy of the speed of light 
through (2.14). 


3 Null test using H(z) and supernovae data 


Given some observational data sets, it is crucial to use a model-independent method to 
reconstruct E(z), D(z), E’(z), D’(z) and D’(z), in order to test the constancy of the speed 
of light following (2.14). Note that here we need also reconstruct the derivatives of the 
functions E(z) and D(z). Here we use the nonparametric approach Gaussian processes [9— 
12] to smooth the data and take the Gaussian processes in Python [12]. The detailed analysis 
and description of the GP method can be found in [12, 13] or in our previous works [14, 15]. 
The analysis of the choice covariance kernel can be also found in [13]. 


3.1 Hubble rate data, Union 2.1 and JLA 


Following [16, 17] we proceed to an analysis based on observational Hubble data compiled 
from several sources, independent of SNela. We combine measurements of H(z) obtained 
with two methods. One is cosmic chronometers, which are mainly passively evolving galaxies. 
There are 21 data points compiled by Moresco et al. [18, 19]. The other is radial baryon 
acoustic oscillations from galaxy clustering in redshift surveys, which give 7 data points of 
Hubble parameters from different experiments [20-23]. We summarize the total 28 data 
points in Table 1. 

To reconstruct D(z), we use two different data, first we adopt the SNela Union 2.1 data 
sets [24], which contain 580 SNela data. Then we employ the most recent SNela catalog 
available: the JLA [25]. We use the binned JLA compilation which shows the same trend as 
using the full catalog itself. For Union 2.1 and the JLA, We have included the covariance 
matrix with both the statistical and systematic uncertainties. We transform the distance 
modulus m — M given in the data set to D using 


m— M +5log 2] — 25 = 5log |(1 + z)D]. (3.1) 
0 
We apply the reconstructions of E(z), D(z) and their derivatives to the null test in Eq. (2.14). 
The result is shown in Fig. 1. 
As expected, from Fig. 1 we can see that for both of the Union 2.1 and the JLA, 
T(z) = 0 lies within lo C.L. of the reconstructed region, which indicates that using current 
observational data there is no signal of variation of c(z). 


3.2 Mock data 


In order to evaluate the constraint on the null test from future observational experiments 
and its ability to probe the small variation of the speed of light from a constant, we create 
mock data to test the ability. Following [5, 6], we choose a general theoretically-motivated 
expression given in [4] as a 


c(a) x co(1 + aa (3:2) 
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69+ 12 18 
83 +8 18 
(Ose 18 
(55 18 
79.69 + 2.32 20 
77+14 18 
82.1 + 4.9 21 
83 + 14 18 
95+ 17 18 


186.5 + 50.4 19 


Table 1. H(z) measurements from different surveys using passively evolving galaxies and radial 


Index z 

1 0.090 

2 0.170 

3 0.179 

4 0.199 

5 0.240 

6 0.270 

7 0.350 

8 0.352 

9 0.400 

10 0.430 

11 0.440 

12 0.480 

13 0.570 

14 0.593 

15 0.600 

16 0.680 

17 0.730 

18 0.781 

19 0.875 

20 0.880 

21 0.900 

22 1.037 

23 1.300 

24 1.363 

25 1.430 

26 1.530 

27 1.750 

28 1.965 

BAO. 

0.4 
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Figure 1. Reconstruction of T(z) from CC+BAO with Union 2.1 (Left) and JLA (Right). The 
shaded blue regions are the 68% and 95% C.L. of the reconstruction. The red line corresponds to co. 
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Figure 2. Reconstructions of T(z) for different models. From top to bottom, Qg = 0,0.16, —0.16, 
respectively. And from left to right, Ac = 0, Ac/co ~ 1%, Ac/co ~ 2%, respectively. The shaded 
blue regions are the 68% and 95% C.L. for the reconstruction. 
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Figure 3. Reconstructions of T(z) for Ac/co ~ 0.1%. The error is one-tenth of the original mock data. 
The shaded blue regions are the 68% and 95% C.L. of the reconstruction. The red line corresponds 
to Co. 


where a = 1/(1 + z) is the scale factor. Note that here we should add a normalization factor 
to make c(a = 1) = cg. Since our model-independent method can probe the constancy of the 
speed of light at every redshift, we can focus on the redshift range where the high quality of 
the data sets lie, i.e. z € [0.2, 0.6]. 

We firstly simulate a data set of 128 points for E(z). Adopting the methodology in [26], 
we use the errors of current observational data. We draw the error from a Gaussian dis- 
tribution: og ~ N(G,¢) with ¢ = (o4 + o_)/2 and e = (0; — o_)/4, where o} and o 
are the two straight lines bound the uncertainties o(z) of the observational E(z) data from 
above and below, respectively. Then E(z)sim is sampled from the Gaussian distribution 
E(Z)sim ~ N(E(2) fia, 72), Where E(z) fig is the theoretical value of the fiducial model. 

For simulated D(z) data, we create mock data sets of future SNela according to the 
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Dark Energy Survey (DES) [27]. The DES is expected to obtain high quality light curves 
for about 4000 SNelas from z = 0.05 to z = 1.2. From Table 14 in [27] we can calculate the 
errors in D and øp and the corresponding numbers of SNela for each redshift bin. At every 
redshift point z, D(z) sim is sampled from the normal distribution D(z) sim ~ N (D(z) fia, oD). 

In our analysis, we produce mock data based on three different cosmological models: 
(1) Ac = 0, the baseline ACDM model; (2) a, = 0.05, n = —0.04, Ac/co ~ 0.5% — 1.5% at 
z € [0.2,0.6]; (3) ac = 0.05, n = —0.09, Ac/co ~ 1.5% — 3% at z € [0.2,0.6]. Additionally, in 
these three models we set Qg = 0, +0.16 and —0.16, respectively, which is significantly large 
to be detected by current data sets using a model-independent method proposed in [15]. If 
the null-tests with different values of the cosmic curvature have the same results, it implies 
that we have indeed dodged the cosmic curvature. Here we emphasize that in the case of 
non-vanishing cosmic curvature, the corresponding Friedmann equation is 

H’ = S, KOU +22. (3.3) 
As illustrated in Fig 2, we can distinguish these three models correctly. Ac/co ~ 1% can be 
detected at ~ 1.50 C.L. and Ac/co ~ 2% or lager can be detected at ~ 30 C.L. The results 
with different values of the cosmic curvature are consistent, which indicates that we indeed 
dodged the cosmic curvature as discussed in Sec. 2. 

Comparing our results to those obtained in the papers [5, 6], we find that DES can not 
provide better improvement in detecting variation of the speed of light for the Ac/co ~ 1% 
case. This is because dodging the cosmic curvature introduces the second derivatives which 
will lead to larger errors on the results. Though we have this disadvantage, our method has 
almost the same ability to detect the Ac/co ~ 1% deviation. If we want to detect smaller 
deviation from cp, we have to use the data sets with higher quality, i.e., bigger number and 
smaller errors of data. Following [6] we also produce the cosmic chronometers H(z) and 
SNela data from DES, but with errors ten times smaller than the expected DES ones. The 
result is shown in Fig. 3. We can see that reducing the errors to one-tenth of the expected 
DES ones will make it possible to detect Ac/co ~ 0.1% at 20 C.L. 


4 Conclusions and discussions 


In this paper, we develop a new method to probe the constancy of the speed of light. By 
dodging the cosmic curvature, we can directly test the speed of light without assuming 
any cosmological model, which should be more natural than the method proposed in [5]. 
Furthermore, we can test c(z) at every redshift covered by the data sets, so that we can 
use real observational data and it is not limited in the constrained region of the redshifts. 
We use H(z) data from the cosmic chronometers and BAO, and then combine them with 
SNIa Union 2.1 and the most recent JLA to give a null test of the speed of light. The result 
indicates that there is no signal of deviation from cg. What we are more concerned about is 
how small variation of the speed of light we can detect using higher quality data sets in the 
future. For this purpose, we create mock data sets based on three fiducial models. We find 
that using the simulated data sets whose errors are obtained from cosmic chronometers and 
DES, we can detect Ac/co ~ 1% at ~ 1.50 C.L. and Ac/co ~ 2% at ~ 30 C.L. If we improve 
the quality of the data, i.e, set the errors to be one-tenth of the expected DES ones, we can 
easily detect a Ac/co ~ 0.1% variation at 20 C.L. The luminosity distance containing more 
a factor ¢(z) outside the integral makes it can detect more information of the speed of light. 
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Dodging the cosmic curvature is crucially important when we want to detect a smaller 
variation of the speed of light because even a very small cosmic curvature can influence the 
results due to the strong degeneracy between them. Our new method developed in this 
paper overcomes such a problem. Although our method has no advantage in the accuracy 
compared to the one in [5] using current data sets, it provides a model-independent way to 
detect variation of the speed of light with more precision when the quality of the data sets can 
be improved in the future. Finally we mention that according to (2.10), in principle, we can 
reconstruct the speed of light, once one knows A(z), B(z), and M(z) from the observational 
data. 

In fact, changing the speed of light at a redshift z leads to a variation in the fine structure 
constant a = e?/(hc), where e is the electron charge and A the reduced Planck constant [28]. 
If the parameters e and A involved in its definition are assumed to be constant, it is easy 
to get Ac/co = —Aa/a, which impacts the measurement of the redshift z of the object. 
Recent analysis in Ref. [29] implies that the variation in the fine structure constant is very 
small, at least Aa/a < 1074. Therefore, we can ignore the effect of the variation in a on the 
measurement of the redshift. In principle, a large variation in the speed of light can still be 
compatible with such orders of magnitudes if the other parameters are allowed to vary. 

It is known that the distance-duality between the luminosity distance and the angular- 
diameter distance is violated by non-conservation of photon number [30, 31]. Of course, such 
a violation of the distance-duality results in the deviation of T(z) in Eq. (2.14) from zero. 
We have to emphasize that only the change of the speed of light has been considered in this 


paper. 
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